Surface decorated silicon nanowires: a route to high-ZT thermoelectrics 



o 
o 



Trocls Markussen, 1 Antti-Pekka Jauho, 1,2 and Mads Brandbyge 1 

1 Department of Micro- and Nanotechnology, Technical University of Denmark, 
DTU Nanotech, Building 345 East, DK-2800 Kgs. Lyngby, Denmark 
1 Department of Applied Physics, Helsinki University of Technology, P.O. Box 1100, FIN-02015 HUT, Finland 

(Dated: July 13, 2009) 

Based on atomistic calculations of electron and phonon transport, we propose to use surface 
decorated Silicon nanowires (SiNWs) for thermoelectric applications. Two examples of surface 
decorations are studied to illustrate the underlying ideas: Nanotrees and alkyl functionalized SiNWs. 
For both systems we find, (i) that the phonon conductance is significantly reduced compared to the 
electronic conductance leading to high thermoelectric figure of merit, ZT, and (ii) for ultra-thin 
wires surface decoration leads to significantly better performance than surface disorder. 
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Recent ground-breaking experiments indicate that 
rough silicon nanowires (SiNWs) can be efficient ther- 
moelectric materials although bulk silicon is not [l|, 
they conduct charge well but have a low heat conductiv- 
ity. A measure for the performance is given by the figure 
of merit ZT = G e S 2 T/n, where G e , S, and T, are the 
electrical conductance, Secbeck coefficient, and (average) 
temperature, respectively. The heat conduction has both 
electronic and phononic contributions, k = K e -\-Kph,. Ma- 
terials with ZT ~ 1 are regarded as good thermoelectrics, 
but ZT > 3 is required to compete with conventional re- 
frigerators or generators [|[. Recent theoretical works 
predict ZT > 3 in ultra-thin SiNWs @, H, @. The high 
performance SiNWs in Ref. [l[ were deliberately pro- 
duced with a very rough surface, and the high ZT is 
attributed to increased phonon-surface scattering which 
decreases the phonon heat conductivity, while the elec- 
trons are less affected by the surface roughness. The ex- 
traordinary low thermal conductivity measured in rough 
SiNWs is supported be recent calculations @, [1] ■ Surface 
disorder will, however, begin to affect the electronic con- 
ductance significantly in very thin wires Q and thereby 
reduce ZT jl, [(| . 

For defect-free wires, the room temperature phononic 
conductance scales with the cross-sectional area (lf3 |. 
Contrary, for ultra-thin wires with diameters D < 5 nm, 
the electronic conductance is not proportional to the 
area: it is quantized and given by the number of va- 
lence/conduction band states at the band edges. De- 
creasing the diameter to this range thus decreases n p h 
while keeping G e almost constant. An optimally designed 
thermoelectric material would scatter phonons but leave 
the electronic conductance unaffected, even for the small- 
est wires. Instead of introducing surface disorder, which 
affects the electronic conductance in the thin wires, it 
is interesting to speculate whether other surface designs 
could lead to improved thermoelectric performance. Lee 
et al. [TT| proposed nanoporous Si as an efficient thermo- 
electric matrial. Blase and Fernandez-Serra [l2j recently 
showed that covalent functionalization of SiNW surfaces 
with alkyl molecules leaves the electronic conductance 
unchanged, because the molecular states arc well sepa- 



rated in energy from the nanowire bandedges, see Fig. 
[T] If the alkyl molecules scatter the phonons, such func- 
tionalized SiNWs would be candidates for thermoelec- 
tric applications. Experimental alkyl functionalization 
of SiNWs was reported in Ref. [l3| . 




Alkyl molecule 



Main wire 




Energy 



Alkyl 

functionalized 
SiNW 



<110> 





Transport direction 



Nanotree 




Side view 



End view 



FIG. 1: Sketch of electronic density of states (top) in the alkyl 
functionalized SiNW (middle) and in the nanotree (bottom). 
The trunk in the nanotree is broader than the branch and thus 
has a smaller bandgap. Also the alkyl molecular states are 
located deep inside the bands. For both structures electrons 
and holes close to the band edges are therefore only weakly 
scattered while phonons are strongly scattered. 

Another possible surface design are branched SiNWs, 
so-called nanotrees. Nanotrees have been synthesized in 
III- VI semiconductors [H,[IH and in silicon (TEE?!- Tne 
stability and electronic structure of silicon nanotrees have 
recently been addressed theoretically [H, Eil ■ The thin- 
ner branches will have a larger band gap than the main 
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wire ("trunk"), see Fig. [TJ Close to the band edges, the 
electronic scattering is therefore weak. 

The presence of an alkyl molecule or a nanowire branch 
leads both to a reduction, Are, of the thermal conduc- 
tance and a reduction, AG, of the electronic conduc- 
tance. In this letter we show that at room temperature 
(RT) the ratio Are/AG > 50 for the alkyl functionalized 
SiNWs, and Are/ AG > 20 for a nanotree. By engineer- 
ing the SiNW surfaces it is thus possible to reduce the 
phonon conductance while keeping the electronic conduc- 
tance almost unaffected. We thus propose such surface 
decorated SiNWs as promising candidates for nanoscale 
thermoelectric applications. 

Systems. We consider two specific systems shown in 
Fig. [1] The first is an alkyl functionalized SiNW with a 
wire diameter of 12 A, and with the wire oriented along 
the (110) direction. The alkyl (C„H 2rl +i) is attached 
to the H-passivated nanowire replacing a H atom. The 
second system is a nanotree, where a small diameter (12 
A) branch is attached to a larger diameter (20 A) trunk. 
The trunk is oriented in the (110) direction while the 
branch is oriented along the (110) direction, and is thus 
perpendicular to the trunk. The length of the branch, 
Lb, is varied. 

Methods. The electronic Hamiltonian, H, and overlap 
matrix, S, of the alkyl functionalized SiNWs are obtained 
from local orbital DFT calculations [13, [H|. The calcula- 
tions are performed on super-cells containing 7 wire unit 
cells with the alkyl molecule bound to the middlemost 
unit cell, as shown in Fig. [T] 

For the nanotrees, we use a tight-binding (TB) model 
since these systems contain > 1100 atoms, too many for 
our DFT implementation. The electronic TB Hamilto- 
nian describing the nanotree is calculated using a 10 band 
sp 3 d 5 s* nearest-neighbor orthogonal TB parametriza- 
tion [2^, HH . We recently applied the same TB model 
to study thermoelectric properties of surface disordered 
SiNWs @]. 

The phononic system, characterized by the force con- 
stant matrix, K, is described using the Tersoff empirical 
potential (TEP) model [H,[25[ for both the nanotree and 
the functionalized SiNW. For pristine wires, we have re- 
cently shown that the TEP model agrees well with more 
elaborate DFT calculations [lj| ■ We limit our description 
to the harmonic approximation, thus neglecting phonon- 
phonon scattering. The harmonic approximation is al- 
ways valid at low temperatures. In bulk Si, the room tem- 
perature anharmonic phonon-phonon relaxation length 
at the highest frequencies is \ a {^>max) ~ 20 nm and in- 
creases as A a oc lo~ 2 at lower frequencies [26j |. Experi- 
mental studies of silicon films [I?! showed that the effec- 
tive mean free path of the dominant phonons at room 
temperature is ~300 nm. For relatively short wires with 
lengths L < 100 nm the anharmonic effects thus seem to 
be of limited importance, and the harmonic approxima- 
tion is expected to be good. 

We calculate the electronic conductance from the elec- 
tronic transmission function, %Ae). This is obtained from 



the H, S-matriccs following the standard non-equilibrium 
Green's function (NEGF)/Landauer setup, where the 
scattering region (i.e. the regions shown in Fig. [IJ is 
coupled to semi- infinite, perfect wires [28|. The elec- 
tronic quantities in the ZT formula can be written 
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Here f(e,fi) = 1/ (exp [(e — fi)/kgT] + 1) is the Fermi- 
Dirac distribution function at the chemical potential pL. 

The phonon transmission function, Tpi^u), at fre- 
quency u) is calculated in a similar way as the elec- 
tronic transmission with the substitutions H — ► K and 
eS — > lj 2 M, where M is a diagonal matrix with the 
[33l Hi)]. The phonon thermal con- 
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FIG. 2: Hole and electron transmissions for pentyl function- 
alized, D — 12 A SiNW (a)-(b) and for a nanotree with 
D = 20 A trunk and D = 12 A branch (e)-(f). Panels (c)- 
(d) and (g)-(h) show the PDOS on the wire and on the pentyl 
and branch, respectively. The energy scales are relative to 
the valence band edge, E v for the holes (left column) and 
relative to the conduction band edge, E c for electrons (right 
column). The DFT band gap of the D = 12 A wire is 1.65 eV, 
while the TB band gap of the D = 20 A wire is 1.77 eV. 

Charge transport. Figures [2] (a) and (b) show the 
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calculated hole and electron transmissions for a pentyl 
(CsHxi) functionalized SiNW. Notably, the transmission 
is nearly perfect close to the bandedges. The average re- 
duction of the transmission in the first hole and electron 
conductance plateaus are 2% and 0. 2%, respectively, in 
agreement with the findings of Ref. [12j ]. Figures 2 (c) 
and (d) show the projected density of states (PDOS) on 
the wire and on the pentyl. The high transmission re- 
gions in panel (a) and (b) are seen to correspond with 
regions of vanishing PDOS on the pentyl molecule. Like- 
wise, at energies in the valence band where scattering is 
observed, there is a relatively large PDOS at the pentyl. 

Figures 2 (e) and (f) show the transmission through 
the nanotree. The branch length is 15.4 A. Again, the 
transmission close to the band edges is nearly perfect, 
with a reduction of 2% and 0.9% for holes and electrons, 
respectively. Figures 2 (g) and (h) show the PDOS on 
the main wire and on the branch. We again observe 
a correspondence between perfect transmission and low 
PDOS on the branch. 

The almost perfect transmissions close to the band 
edges can be qualitatively understood from the schematic 
drawing in Fig. Q] (top). The HOMO and LUMO level 
of the pentyl are located deep inside the bands [T3] and 
the molecular states are thus not accessible for electrons 
or holes close to the band edges. For the nanotree, the 
branch has a smaller diameter and thus a larger bandgap. 
Electrons or holes in the trunk, with energies close to the 
band edges, are not energetically allowed in the branch 
and therefore do not 'see' the branch. In addition to 
the energy considerations, the spatial distribution of the 
Bloch state also plays a role: the first valence and con- 
duction band Bloch states of the main wire have more 
weight in the center of the wire than at the edge @ . 

Thermal transport. Figure [3] shows the temperature 
dependence of the thermal conductance ratios k/kq, 
where kq is the pristine wire thermal conductance, which 
in the low energy limit equals the universal thermal con- 
ductance quantum, kq(T) = 4(tt 2 k B T / 3h) [35j. Panel 
(a) shows the ratios for wires with alkyles, C„H2„+i, with 
different lengths, n = 3, 5, 7. The thermal conductance 
at RT is reduced by ~ 10%, and the overall behavior 
does not depend on the alkyl length. The inset shows 
the phonon transmission at low phonon energies. Note 
the resonant dips in the transmission, where exactly one 
channel is closed yielding a transmission of three. These 
dips are associated with an increased local phonon den- 
sity of states at the alkyl molecule at the resonant ener- 
gies, corresponding to a localized vibrational mode. Such 
Fano-like resonant scattering is well-known from elec- 
tron transport jUJ . A phonon eigenchannel analysis (37| 
shows that the transmission dips are due to a complete 
blocking of the rotational mode in the wire. The corre- 
sponding localized alkyl phonon mode is a vibration in 
the plane perpendicular to the wire axis. 

Panel (b) in Fig. [3] shows the thermal conductance ra- 
tio for nanotrees with different branch length, Lb- There 
is only a weak dependence on Lb at low temperatures, 
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FIG. 3: Thermal conductance ratio k/ko for (a) pentyl func- 
tionalized SiNWs with different diameters, and (b) for nan- 
otrees with different branch lengths, Lb- The nanotree trunks 
are again D — 20 A with 12 A diameter branches. The in- 
sets show the phonon transmission function at low energies. 
Fano-like resonant scattering is observed in both systems. 



and at RT the four curves basically coincide showing a 
thermal conductance reduction of 17% of the nanotree 
compared to the pristine wire. Again we observe reso- 
nant transmission dips for the nanotrees. Two channels 
- the rotational and one flexural mode - close completely 
at the resonance due to two quasi-localized vibrational 
modes in the branch. These phonon backscattering res- 
onances are responsible for the dip in the k/kq ratio 
around T = 10 K. Notice that all the conductance ra- 
tios approach unity in the low temperature limit. This 
is because the four acoustic modes transmit perfectly in 
the limit uj -» [HI]. 

We may vary the thermoelectric figure of merit, ZT, 
by varying the chemical potential. Typically ZT displays 
a maximum for /1 close to the band edge 0, El . Figure 
2] shows the maximum ZT values for the pentyl func- 
tionalized SiNW (squares), the nanotree (circles), and 
surface disordered SiNWs (triangles), where disorder is 
modeled by introducing surface silicon vacancies. The 
diameter of the surface disordered wire is D = 20 A and 
it is oriented along the (110) direction. The calculational 
details are given in Ref. Q. The curves show ZT as a 
function of the number (A) of pentyl molccules/nanotrcc 
branches/silicon vacancies. In calculating ZT vs N, 
we have assumed that the transmission, T/v ; through a 
longer wire with e.g. N pentyl molecules randomly cov- 
ering the surface can be obtained from the single-pentyl 
transmission, Ti as T^ 1 = T^ 1 + A(T X _1 — T^ 1 ), where 
7o is the pristine wire transmission. The term in paren- 
thesis corresponds to a scattering resistance of a single 
pentyl molecules. This averaging method has recently 
been validated in the quasi-ballistic and diffusive regimes 
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for both electron and phonon transport [1, [H, [3l| . 
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FIG. 4: Thermoelectric figure of merit, ZT, for p-type (a) 
and n-type (b) wires. N is the number of pentyl molecules 
(squares), nanotree branches (circles) and silicon surface va- 
cancies (triangles) in the wire. 

Figure |4] shows that increasing the number of scatter- 
ing centers, i.e. the number of pentyl molecules or nan- 
otree branches increases the ZT for both hole transport 
(a) and electron transport (b) . In the case of holes in the 
pentyl functionalizcd SiNWs, the ZT reaches an almost 
constant level of ZT = 0.7 at N = 40, but in all other 
cases, ZT increases throughout the range. Increasing the 



density of moleculcs/nanotrec branches or increasing the 
length of the wire will thus increase the thermoelectric 
performance. The reason is that the electrons (holes) 
are less affected by the surface modifications than the 
phonons, as also seen in Figs. [2] and [3] The surface dis- 
ordered wires (triangles) show an increasing ZT vs N but 
at values significantly lower than the two other surface 
modified wires. 

Discussion. A number of idealizations have been 
made in our calculations, and we next assess their sig- 
nificance. The structures we have considered represent 
plausible choices, dictated by computational limitations, 
but do not necessarily match quantitatively real struc- 
tures. Thus, for example, a surface decorated SiNW will 
also be rough, and one should consider the combined ef- 
fect of all scattering mechanisms. We have not carried 
out optimizations neither with respect to the attached 
molecules nor with respect to the geometry of the nan- 
otrees. Electron-phonon and phonon-phonon scattering 
will affect both the electronic and thermal conductances 
and the obtained ZT values We do not expect to 
reach quantitative agreement with experiment but be- 
lieve to have identified important trends: In SiNW based 
thermoelectrics, surface decorations in terms of added 
molecules or nanowire branches seem to be a better ap- 
proach than surface disorder in the ultra-thin limit. 
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